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LONG-TERM  GOALS 

The  long-term  objective  is  to  contribute  to  the  development  of  the  components  of  limited  area,  open 
boundary,  coastal  nowcast/forecast  systems  that  will  resolve  the  time  and  length  scales  of  the  relevant 
physical-biological  dynamics  in  shallow  coastal  environments. 

OBJECTIVES 

Our  objectives  are  to  develop  the  methodology  for  bioluminescence  potential  (BL)  and 
bioluminescence  leaving  radiance  (BLw)  predictions  on  scales  to  1-5  days,  and  to  understand  the 
coupled  bio-optical  and  physical  processes  in  the  coastal  zone  that  governs  the  variability  and 
predictability  of  bioluminescence. 

APPROACH 

Approach  is  based  on  joint  studies  of  the  marine  bioluminescence  potential  (BL)  and  Inherent  Optical 
Properties  (IOPs)  over  relevant  time  and  space  scales.  The  physical  model  is  based  on  the  Navy 
Coastal  Ocean  Model  (NCOM),  the  biochemical  model  simulates  dynamics  of  two  sizes  of 
phytoplankton,  zooplankton,  nutrients,  ammonium,  and  detritus  pools  (Shulman  et  ah,  2011).  The  BL 
potential  model  is  based  on  advection-diffusion-source  model,  with  velocities  and  diffusivities  taken 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2012 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Dynamic  Modeling  of  Marine  Bioluminescence  and  Night  Time  Leaving 
Radiance 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory  Stennis  Space  Center,  MS  39529 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

SAR 

18.  NUMBER 
OF  PAGES 

5 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


from  the  physical  model  (Chai  et  ah,  2002;  Fujii  et  ah,  2007;  Shulman  et  ah,  2011).  The  proposed 
research  is  being  significantly  leveraged  by  the  interdisciplinary  and  multi-institutional  modeling  and 
field  efforts  of  the  NRL  BIOSPACE  and  MURI  ESPRESSO  projects.  Bio-optical,  physical 
observations  from  the  following  field  programs  are  being  used  in  this  study:  AOSN  I  and  II  (Moline, 
2007  and  2009,  Shulman  et  ah,  2003,  2005,  2009);  bioluminescence  observations  from  Monterey  Bay 
surveys  conducted  by  Dr.  Haddock  (Augusts  of  2000,  2002,  Decembers  of  2002  and  2003  and  March 
of  2004);  NRL  BIOSPACE-ESPRESSO  May- June  of  2008  experiment;  NRL  BIOSPACE/MB ARI 
CANON  experiment  in  October  2010.  We  design  predictive  experiments  when  the  model  is  initialized 
by  using  a  set  of  BL  potential  observations.  After  that  we  conduct  forward  in  time  simulations  without 
assimilation  of  new  observations,  in  this  case  the  experiments  simulate  the  forecasting  of  the  BL 
potential.  The  forecasts  are  compared  to  the  non-assimilated  observations.  Bio-optical,  physical 
observations  are  being  used  also  for  verification  and  interpretation  of  our  numerical  experiments 
results. 

WORK  COMPLETED 

Dynamical,  predictive  biochemical,  physical  and  bioluminescence  potential  (BL)  models  are  combined 
into  a  methodology  for  short-term  modeling  and  forecasting  of  BL  potential,  Inherent  Optical 
Properties  (IOPs)  and  the  nighttime  water-leaving  radiance  (BLw)  due  to  stimulated  bioluminescence 
at  depth. 

We  demonstrated  that  observed  bioluminescent  dino flagellates’  avoidance  of  advection  during  the 
upwelling  development  can  be  explained  by  their  vertical  swimming  ability.  This  mechanism  also 
helps  explain  the  general  persistence  of  dino  flagellates  in  the  upwelling  shadow  area.  We  published 
JGR  Oceans  paper  on  this  subject.  The  paper  was  selected  by  AGU  Editors  as  AGU  “Research 
Spotlight”. 

Project  results  were  presented  at  the  ONR  Marine  Mammals  &  Biology  Program  Review,  April  16-18, 
2012,  Washington,  DC. 

Book  Chapter  (which  includes  results  from  this  project)  has  been  submitted. 

RESULTS 

In  the  upwelling-driven  system  of  Monterey  Bay,  C  A,  results  showed  that  the  offshore  water  masses 
with  the  subsurface  layer  of  bioluminescent  zooplankton  were  replaced  by  water  masses  advected  from 
the  northern  coast  of  the  bay  with  a  relatively  high  presence  of  mostly  non-bioluminescent 
phytoplankton  (Shulman  et  ah,  2011).  Offshore  observations  show  a  deeper  BL  potential  maximum 
below  the  surface  layers  of  high  chlorophyll  and  backscatter  values  during  the  earlier  stages  of 
upwelling  development.  Later,  the  observed  deep  offshore  BL  potential  maximum  disappeared  and 
became  a  shallower  and  much  weaker  signal.  Our  analysis  showed  that  during  the  upwelling, 
bioluminescent  dinoflagellates  from  the  northern  part  of  the  Bay  were  able  to  avoid  advection  by 
strong  southward  currents  developed  during  the  upwelling  event.  Predictive  physical,  biochemical  and 
bioluminescence  intensity  models  were  used  to  interpret  the  observed  dynamics  during  the  upwelling. 


2 


Figure  1.  (right)  Water-leaving  radiance  at  the  surface  due  to  stimulation  of  the  modeled  BL 
potential  at  5m  depth.  The  modeled  BL  potential  at  5m  depth  (middle  panel),  and  a  sum  of  a 
(absorption)  and  bh  (backscattering)  averaged  from  the  depth  of  BL  potential  stimulation  (5m) 

to  the  surface  (left  panel). 

The  results  (Figure  1)  show  high  values  of  BL  potential  and  estimated  BLw  in  areas  along  the 
coastline,  where  high  values  of  BL  potential  were  observed  and  model-predicted.  At  the  same  time,  the 
model  showed  high  values  of  BL  potential  along  the  entrance  to  the  Bay  (Figure  1,  Shulman  et  al., 
2011),  which  was  not  in  agreement  with  observations.  In  the  model,  the  BL  potential  dynamics  were 
controlled  by  advective  and  diffusive  processes  only,  and  as  it  was  speculated  in  Shulman  et  al,  2011, 
the  lack  of  modeling  of  behavioral  dynamics  of  bioluminescent  organisms  (as  for  example,  vertical 
migration  of  dino flagellates  to  avoid  advection  losses)  is  responsible  for  the  model’s  inability  to 
predict  the  observed  weakening  of  the  BL  potential  along  the  entrance  to  the  Bay.  In  Shulman  et  al., 
2012,  we  tested  the  hypothesis  that  vertical  swimming  behavior  explains  observed  ability  of 
dino  flagellates'  to  avoid  advection  by  strong  currents.  The  dynamics  of  dino  flagellates  was  modeled 
with  the  tracer  model  where  the  dinoflagellate  population  is  modeled  as  a  concentration,  and  vertical 
swimming  velocity  is  introduced  into  the  tracer  advective-diffusive-source  model.  Three  swimming 
behaviors  were  considered:  sinking,  swimming  to  the  target  depth  and  diel  vertical  migration. 
Swimming  velocities  in  all  cases  were  considered  in  the  range  of  documented  velocities  for  the 
dinoflagellates  species  observed  during  the  upwelling  development  in  the  Monterey  Bay.  Results 
demonstrated  that  through  swimming  behavior,  dinoflagellates  avoid  complete  advection  out  of  the 
Bay  during  upwelling  events  (Shulman  et  al.,  2012).  With  a  modeled  swimming  velocity  of  20m/day 
(a  reasonable  estimate  at  half  the  observed  maximum)  40%  of  the  dinoflagellates  population  was 
advected  from  the  northern  part  of  the  Bay  compared  to  no  swimming.  This  is  in  agreement  with  the 
observed  mean  BL  potential  ratio  of  0.45  at  the  Bay  entrance  compared  to  the  northern  part  of  the  Bay. 
While  some  of  the  salient  features  of  short-tenn  changes  in  BL  potential  can  be  predicted  and 
explained  with  the  modeling  of  advective-diffusive  processes,  it  was  demonstrated  that  the  modeling  of 
BL  sources  and  sinks  is  needed  to  reproduce  the  observed  spatial  and  temporal  variability  of  the  BL 
even  on  short-time  scales.  It  is  especially  valid  in  situations  when  swimming  behavior  of 
bioluminescent  plankton  impacts  the  BL  potential  distribution.  Continued  modeling  efforts  such  as 
these  will  gauge  the  measure  of  bioluminescence  as  a  tool  for  integrating  ecosystem  infonnation, 
evaluate  dynamical  optical  properties  in  the  ocean,  and  help  in  short-term  prediction  of  oceanographic 
conditions,  including  BL  potential  and  BLw. 
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IMPACT/APPLICATIONS 


Prediction  of  the  location,  timing  and  intensity  of  bioluminescence  potential  is  critical  for  numerous 
naval  operations  including  preventing  detection  of  covert  operations  involving  submarines,  Swimmer 
Delivery  Vehicles  and  AUVs,  as  well  as  in  aiding  detection  of  enemy  incursions.  At  present,  the  Navy 
does  not  have  capability  to  forecast  BP  potential  and  night  time  water  leaving  radiance.  The  proposed 
research  aims  to  develop  a  methodology  for  bioluminescence  potential  and  bioluminescence  leaving 
radiance  predictions  on  scales  to  1-5  days. 

TRANSITIONS 

None 

RELATED  PROJECTS 

NRL,  RO  "  Bio-Optical  Studies  of  Predictability  and  Assimilation  for  the  Coastal  Environment 
(BIOSPACE)”  (PI:  I.  Shulman) 

I.  Shulman  is  PI  of  the  NRL  project  with  objectives  to  improve  understanding  of  the  variability  and 
predictability  of  the  underwater  light  and  bio-optical,  physical  properties  on  time  scales  of  1  to  5  days. 
NRL  coupled  models  and  predictions  of  physical  bio-optical  properties  (including  IOPs  and  BP)  are 
used  in  our  project. 

The  Multidisciplinary  University  Research  Initiative  (MURI)  project  “Rapid  Environmental 
Assessment  Using  an  Integrated  Coastal  Ocean  Observation-Modeling  System  (ESPRESSO)”  (Pis:  O. 
Schofield,  S.  Glenn,  J. Wilkin,  G.  Gawarkiewicz,  R.  He,  D.  McGillicuddy,  K.  Fennel,  M.  Moline). 

Objectives  of  the  MURI  project  are  focused  on  the  development  of  a  data  assimilative  physical-optical 
modeling-observation  system  capable  of  improving  predictive  skill  for  forecasting  ocean  color  and 
improving  physical  models  by  using  ocean  color.  M.  Moline  is  Co-PI  of  the  project,  and  NRL 
BIOSPACE  and  MURI  project  have  similar  objectives  and  there  are  ongoing  collaborations  between 
projects. 
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